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We address the question of how one can modify the inflationary tensor spectrum without changing 
at all the successful predictions on the curvature perturbation. We show that this is indeed possible, 
and determine the two quadratic curvature corrections that are free from instabilities and affect only 
the tensor sector at the level of linear cosmological perturbations. Both of the two corrections can 
reduce the tensor amplitude, though one of them generates large non-Gausslanity of the curvature 
perturbation. It turns out that the other one corresponds to so-called Lorentz-violating Weyl gravity. 

In this latter case one can obtain as small as 65% of the standard tensor amplitude. Utilizing this 
effect we demonstrate that even power-law inflation can be within the 2a contour of the Planck 
results. 

PACS numbers: 98.80.Cq 


I. INTRODUCTION 

Inflation [1, 2] plays a crucial role in cosmology of 
the very early Universe. In particular, single-field slow- 
roll models of inflation generically produce nearly scale- 
invariant, adiabatic, and Gaussian curvature perturba¬ 
tions as the seeds for cosmic structure [3]. The theo¬ 
retical prediction matches observational results e.g. of 
the Planck experiments fairly well [4—7]. During infla¬ 
tion primordial tensor modes (gravitational waves) are 
generated as well. The tensor amplitude is convention¬ 
ally parametrized by the tensor-to-scalar ratio, r, and 
the Planck constraint on r is given by r < 0.10 (95% 
C.L.) [7]. Some of the single-field slow-roll models pre¬ 
dict larger tensor modes, and hence have been excluded 
by current observations. One would then ask whether 
one can reduce the tensor amplitude somehow to save 
such models. This is the question which we discuss in 
this paper. 

General relativity is an underlying assumption of stan¬ 
dard inflation models, and nonstandard dynamics of the 
tensor modes can be obtained by modifying this gravita¬ 
tional sector. In doing so, one generically expects that 
the behavior of the scalar perturbations is also modified. 
This is however what we want to avoid, because the stan¬ 
dard inflationary predictions on the scalar perturbations 
are so successful. In this paper, we therefore explore the 
possibility of modifying only the tensor modes and try 
to retain the same structure of the scalar sector as in 
general relativity, in order not to spoil the remarkable 
agreement between the standard theoretical predictions 
of the scalar perturbations and observations. 

It is natural to consider quadratic curvature terms in 
the action beyond general relativity since such correc¬ 
tions are expected to arise as signatures of new physics 
at high energies. Below we look for quadratic curvature 
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terms that modify only the tensor sector of cosmolog¬ 
ical perturbations without introducing any pathologies 
such as ghost instabilities. It turns out that there are 
two independent combinations of the curvature tensors 
fulfilling the above requirements. Both combinations do 
not change the quadratic action for the scalar perturba¬ 
tions, and one of them has no impact on the cubic action 
as well. The resultant quadratic curvature terms are not 
of the form TZ^, , etc. which are familiar in the 

literature [8, 9], but they have nontrivial coupling to the 
derivative of the inflaton field. We study in detail how 
the tensor amplitude and tilt are modified, and discuss 
the implications for observations. 

The organization of this paper is as follows. We deter¬ 
mine the two possible quadratic curvature terms satisfy¬ 
ing our requirement in the next section. In Sec. Ill we 
evaluate the modified amplitude and tilt of the primor¬ 
dial tensor modes, and then present the implications for 
observations. We draw our conclusions in Sec. IV. 


II. CONSTRUCTION OF THE LAGRANGIAN 

The theory we consider is described by the action 

S = S'eH + S^-\- S'higher, (1) 

where S'eH is the Einstein-Hilbert term, 

-S'eH = ^ y (2) 

is the action of the inflaton field, 

S0 = J d'^Xy/^P{(j),d^cj)df,cj)), (3) 

and Shigher represents higher curvature corrections, 
Shigher = IJ P- + . . . ^ (4) 
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The simplest Lagrangian for the inflaton field would be 
of the form P = — V{(j)), but here we do not 

need to specify the concrete form of P. 

It is known that typical higher curvature terms like 
the one presented in Eq. (4) give rise to new propagating 
degrees of freedom which are plagued by (ghost) insta¬ 
bilities [10]. In this paper, we carefully construct higher 
curvature terms so that the resultant theory is free from 
such dangerous degrees of freedom. Among such healthy 
theories we are interested in those in which the dynamics 
of tensor perturbations is modified while the scalar sector 
of cosmological perturbations is left unchanged. To find 
the higher curvature terms having those properties, we 
have to go beyond the familiar curvature invariants such 
as and , and consider the terms ob¬ 

tained by contracting with the unit normal to constant 
(j) hypersurfaces, 

and the induced metric, 

^pu — Qfiu T (6) 

e.g., Tlpupa'R-p.'u'l''"''Y’’'■ This possibility 
was demonstrated in the context of Weyl gravity in 
Ref. [11]. 

Focusing on quadratic curvature corrections, we are go¬ 
ing to identify the terms in the Lagrangian fulfilling the 
above requirements in the following way. The basic idea 
here is along the same line as taken in Refs. [12, 13]. We 
start by performing the Arnowitt-Deser-Misner (ADM) 
decomposition, taking constant (j) hypersurfaces as con¬ 
stant time hypersurfaces, as the dynamics of cosmological 
perturbations is more transparent in the ADM language. 
The metric is written in terms of the ADM variables as 

ds2 = + 7y (dx* -h iVMt) {dx^ + Wdt). (7) 

The possible quadratic curvature terms in the La¬ 
grangian are exhaustively written in terms of the three- 
dimensional geometric quantities as 

X {K\ 

K^R, ■■■ , • • • } , ( 8 ) 

where Kij and Rij are the extrinsic and intrinsic cur¬ 
vature tensors of the constant tp hypersurfaces, respec¬ 
tively, Di stands for the covariant derivative with respect 
to 7 ij, and ellipses are used to indicate analogous terms 
whose indices are contracted in different ways. We dis¬ 
card from the above candidates the terms containing time 
derivatives of the extrinsic curvature, because higher time 
derivatives of the metric imply the appearance of addi¬ 
tional propagating degrees of freedom other than p and 
two tensor modes, signaling instabilities.^ 


^ If one has only K the theory is not necessarily unstable, as is 


Let US consider cosmological perturbations, 

N = 1 + SN, Ni = diX + Xi7 7ij = , 

(9) 

where C is the curvature perturbation on the uniform 
(j) hypersurfaces, hij is the transverse and traceless ten¬ 
sor perturbation, and Xi is the transverse vector pertur¬ 
bation. Let us concentrate on the scalar sector for the 
moment. To first order in perturbations, the extrinsic 
curvature is given by 

K/ = H5i + ^(577(5/ + SK^, (10) 

o 

with 

5K = -?>H5N+ ^d^X: (11) 

5Ki = - ^ ( 9 . 9 ^ - X, (12) 

and the intrinsic curvature is 

5Ri = (13) 

The perturbation of the extrinsic curvature tensor has 
been decomposed into its trace and traceless parts. Using 
those quantities as building blocks, one can construct the 
following two combinations of the form listed in Eq. (8) 
for which the scalar-type variables are canceled out after 
integration by parts, 

2di6Kjkd^SK^^ - SK^k, (14) 

and 

SR^jSR^^ - IsR"^, (15) 
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at quadratic order in perturbations. No other combina¬ 
tions can be found with vanishing scalar-type variables. 
Now including vector and tensor perturbations we have 

2di5kjkd^5K^^ - ddiSK^'^d^SKjk 

= i(®‘U)' + 5^Ux,)b (16) 

SR,j6R^^ - ^6R^ = ^ • (17) 

Both of the two possible quadratic terms for hij with four 
derivatives are obtained, while we successfully exclude 
hfj which would cause Ostrogradski ghosts. Since there 
is no kinetic term for Xi here and in S'eh, the vector 


illustrated by the example of the model. This however adds 
an extra scalar degree of freedom modifying the scalar sector of 
cosmological perturbations. For this reason we avoid any time 
derivatives of Kij. 
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perturbation is not dynamical. We therefore ignore the 
vector sector in this paper. 

Having thus written the quadratic Lagrangian for per¬ 
turbations in terms of the geometric quantities, it is al¬ 
most straightforward to determine the full nonlinear La¬ 
grangian in the ADM form as 

, (18) 


with and hence it is Lorentz violating in the same 
sense that £i is. 

In Cl and £2 in the ADM form one may consider 
time-dependent M. This translates to the ^-dependent 
coupling in the covariant language. However, since (j) is 
slowly rolling, it is reasonable to assume that M is only 
weakly time dependent. For simplicity we treat M as 
constant in the following. 

III. THE TENSOR AMPLITUDE 


where ATy := ATy — (l/3)Ar7y is the traceless part of 
the extrinsic curvature. Note that £( is one of the can¬ 
didates; in fact, we have different choices that reduce 
to Eq. (14) after integration by parts at the level of the 
quadratic Lagrangian for perturbations. Among them we 
adopt 


£1 = 

- 2D,kjkD^k^^) (20) 

rather than £(. What is particular to £1 is that it can 
be written as a square of some tensor as 

Cl = ( 21 ) 

where 

hLijfc = SDfiATjjfc + (22) 

It is clear that the scalar perturbations do not partici¬ 

pate in Wijk at first order. This means that £1 does not 
modify the scalar sector both at quadratic and cubic or¬ 
der. In other words, the prediction for non-Gaussianity of 
the curvature perturbation, as well as that for the power 
spectrum, remains the same in the presence of £ 1 . This 
is however not the case for £( and £ 2 . 

The covariant form of the Lagrangian can be recovered 
by writing the extrinsic curvature as VpUcr 

and making use of the Gauss-Codazzi relations: 


+ K^Ksp - A7A„^,(23) 

= D^K^p-DpK^^. (24) 

We find 


r - r 

Li — , 0| 


M2 


pupa Cp'u'p' 




(25) 


where C^^pa- is the Weyl tensor. Thus, it turns out that 
£1 reproduces the theory studied in Ref. [11]. One can 
repeat the same procedure also for £2 to derive its co¬ 
variant form. However, the covariant expression for £2 
is messy and not so illuminating, so that in this case it 
is better to work in the simpler ADM form. It is worth 
noting that the covariant expression for £2 is also con¬ 
structed by contracting the Riemaim curvature tensor 


In the previous section we have identified the two pos¬ 
sible quadratic curvature terms that make no contribu¬ 
tion to the scalar sector of cosmological perturbations 
at least at linear order. Let us now investigate how the 
amplitude of primordial tensor modes is modified due to 
those terms. For clarity we study each term separately 
below. Actually, we will find that a sizable modification 
from £2 is prohibited because £2 also produces large non- 
Gaussianity. 


A. £1 


First we consider 


‘^'higher — / ^ xCi, 


(26) 


The quadratic action for the tensor perturbations is given 
by [11] 


S = 


8k 


dt(kxk hij - ^{dkhijf + 


(27) 


Each Fourier mode of two polarization states, h^{t) (A = 
-b, x), obeys a second-order evolution equation. We use 
the canonically normalized variable 


- I i:) 


1/2 


^3/2 


2 \ 1/2 


1-b 


4/c^ \ 

JFkk ) 


K, 


(28) 


and omit A when unnecessary. We then have 

/fc+7(t)/fe = 0, (29) 

where 


K' 


uit := --lH^ + 2H 


Y/Y -2H^-H 
l-b4fc2/M2a2 




(1 -b4P/M2a2)2‘ 
We use the WKB solution. 


fk ^ 


1 


y/20Jk 


exp 


uJk{t')dt' 


(30) 


(31) 
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for the short wavelength modes with k^jd} ^ 

In the short wavelength limit, Eq. (30) reduces to 

+ (32) 

From this we see that during inflation the tensor pertur¬ 
bations are stable if 


H <M. 


(33) 


We assume that the background evolution, which is con¬ 
trolled by the inflaton sector 5^, satisfies the condi¬ 
tion (33). Since Eq. (32) gives the estimate 


ujk eH^ 

OjI ^ (M2 - 772)3/2 ’ 


(34) 


where e := —H/H^ <C 1 is the slow-roll parameter, the 
WKB approximation is justihed as long as H is not too 
close to M. 

In the long wavelength limit, k^/a^ <C we 

have 


Ll)2~— -772_^w_ 


(35) 


so that the standard result is recovered on superhorizon 
scales, hk — const. 

Let us compute the power spectrum of the tensor 
modes. 


VT{k) = ^\hk\\ (36) 

In general, Eq. (29) cannot be solved analytically, and 
hence one needs numerical calculations to evaluate the 
power spectrum. However, in the special case of the exact 
de Sitter background one can solve Eq. (29) analytically 
using the hypergeometric functions. This was done in 
Ref. [11], and here we quote their final result: 

2k772 , , ^ 

IPt = iJ/M , 37 

TT^ 


where 


^l(x) 


cosh(7rj^/2) coth(7rzz/2)|r(—1/4 -|- ij^/4)|"*' 

1287r2a:3 


,(38) 


with 1 / := d a;“2 — l. Based on this, one may expect that 
in the case where 77 is varying the power spectrum is 
given by evaluating the de Sitter result (37) at horizon 
crossing. 


PT(fc) 


2k772„ 


(77/M) 


k—aH 


(39) 


as is commonly done in general relativity. 

We numerically solved Eq. (29) in the case of power- 
law inflation, a oc , using the initial condition (31), and 
verified that Eq. (39) reproduces the numerical results 
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FIG. 1: The power spectrum Vt of tensor modes from power- 
law inflation with the correction from C\. The upper line and 
points are for a cx and the lower for a (x , respectively. 
Red points represent the numerical results, while the dashed 
line indicates the analytic estimate (39). The parameters for 
the upper line and points are given by I tend ~ 10 and 

— 0.16, where tend is the time at the end of inflation. 
The parameters for the lower ones are given by = 

10-^ and = 0.74. 



FIG. 2: 72 as a function of H/M. 


very accurately, as shown in Fig. 1. We are thus allowed 
to use the formula (39) for slow-roll inflation. 

The behavior of the function Si(a;) is as follows: Si —>■ 
1 as a; —>■ 0, and Si < 1 for 0 < a: < 0.95. Thus, the 
tensor amplitude is suppressed for 77 < 0.95M. The 
minimum of Si is given by Si ~ 0.65, which occurs at 
X ~ 0.74, and Si diverges as x —?► 1. 

Since the power spectrum of the curvature perturba¬ 
tion remains unchanged in our theory, the tensor-to- 
scalar ratio is given by 


r = 16eSi. (40) 

The tensor tilt, tit '■= dlnT^r/dlnfe, is evaluated as 
2e 


Tlx — 


l-e 


1 + 1. '!>"=■ 


2dln(77/M) 


(41) 


k=aH 
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We see that din Si/din a; < 0 for 0 < x < 0.74 and 
its minimum value is given by dlnSi/dlnx ~ —0.46 at 
X ~ 0.53. This shows that the tensor spectrum is always 
red. 

From Eqs. (40)-(41) it is clear that the consistency 
relation [14] is violated. The deviation from the stan¬ 
dard consistency relation is characterized by the follow¬ 
ing function, 


V := 


1 -I- (l/2)dln.=,i/dlnx 


(42) 


x=H/M 


as —SriT/r ~ 'D\k=aH- In Fig. 2, we plot as a function 
of H/M. We see that the violation depends on the scale 
k, and Fig. 2 tells us its scale dependence. 

Figure 3 illustrates the observational implications of 
the Cl correction by comparing the suppressed tensor 
amplitude with the Planck results. The red stars in the 
figure indicate the case of power-law inflation, assuming 
the maximal suppression (Si = 0.65) at the observed 
scale. Although original power-law inflation (represented 
by the dashed line) is ruled out by observations, it can 
be within the 2a contour with the help of Ci. The same 
applies to other inflation models such as F (x cjP'. Those 
models originally predict large tensor modes, but the Ci 
correction can bring such models to the observationally 
preferred region in the Ug-r plane. 


B. £2 


Next let us consider 

•Shigher = d'‘a;£ 2 - (43) 

Here we added a minus sign so that the tensor perturba¬ 
tions are stable at high momenta. The quadratic action 
for the tensor perturbations is 


may be used for the short wavelength modes, because 
w/'^dwfc/dry <C 1 is always satisfied at large k. 

Only in the case of exact de Sitter inflation for which 
the scale factor is given by a = l/i7(— 77 ), Eq. (45) can be 
solved analytically. The solution that matches Eq. (47) 
at large k is obtained in terms of the Whittaker function 
as [16, 17] 


_ e ^/^^IFi/4rr,3/4(-ia7fc^77^) 
{—2xk‘^T])^^'^ 

where x = H/M. This yields the power spectrum 
Vt = 

TT^ 


(48) 


(49) 


where 


(x) := ^ [e"/(4")x3/2lr(5/4 + i/(4x))l^ 


4 L 


1 -1 


(50) 


In the case of slow-roll inflation, one may evaluate the de 
Sitter result (49) at horizon crossing, k = aH. This can 
also be justified by a numerical calculation. 

One sees that S 2 is a monotonically decreasing func¬ 
tion and S 2 —>■ 1 as X —>■ 0. Therefore, also in this case the 
tensor amplitude is suppressed relative to the standard 
result. Since S 2 oc x“^/^ for x ^ 1, the tensor ampli¬ 
tude could potentially be suppressed to a very small level. 
However, as we see below, this possibility is hindered by 
the generation of large non-Gaussianity. 

In contrast to the case of £ 1 , the cubic action for the 
curvature perturbation is affected by £ 2 . This implies 
that M must be sufficiently large in order to avoid large 
non-Gaussianities in C,. Typically, £2 contains terms such 
as 


£2 


M 2 


adH? 


(51) 


j dtd^xa^ 



(44) 


Using the conformal time defined by dry = dt/a and the 
canonically normalized variable := (4K)“^/^a/i^ in the 
Eourier space, we obtain 


The non-Gaussianity generated by this term is estimated 
to be 

"" Tap ■ 

Requiring that /nl ^ 1, we have 


where 



+ wl{r])vk = 0 , 


:= P 


M 2 a 2 


I d^g 
a dr]'^ 


(45) 


(46) 


This modified dispersion relation has been studied in de¬ 
tail in the literature [15, 16]. 

The WKB solution 


Vk 


y/2.ujk 


exp 



uJkir]')dr]' , 


(47) 


H 

M 


< « 1 . 


(53) 


Therefore, in fact the suppression factor S 2 cannot be 
much smaller than 1. We conclude that the second La- 
grangian £2 does not provide an efficient way of suppress¬ 
ing the tensor amplitude. 

One may consider a combination of the two La- 
grangians, a£i-|- 6 £ 2 - Obviously, this does not change the 
quadratic Lagrangian for the scalar perturbations, and 
to avoid large non-Gaussianities we must require b 1. 
Therefore, to suppress the tensor amplitude most effec¬ 
tively, essentially one can only use £ 1 . 
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FIG. 3: Reduced tensor amplitude and the Planck results [7] in the ris-r plane. Red stars correspond to power-law inflation 
with the jCi correction, assuming the maximal suppression, Hi = 0.65. 


IV. CONCLUSIONS 

In this paper, we have studied inflationary predictions 
of theories with quadratic curvature corrections. We be¬ 
gan by looking for ghost-free quadratic curvature terms 
that retain the same quadratic action for the curvature 
perturbation as in general relativity while modifying the 
dynamics of tensor perturbations. We have shown that 
such curvature terms can indeed be constructed, and de¬ 
termined the two possible combinations (denoted by £i 
and £ 2 ). This was done by using the ADM formalism, 
and recast in a covariant form those corrections contain 
the curvature tensors contracted with the unit normal 
to hypersurfaces on which the inflaton is homogeneous. 
It has turned out that one of the two terms, £ 1 , is in 
fact identical to the one introduced in so-called Lorentz- 
violating Weyl gravity [11]. This term does not change 
the action of the curvature perturbation even at cubic or¬ 
der. The other term, £21 in contrast, modifies the scalar 
sector at cubic order. 

We have investigated the tensor amplitude in the pres¬ 
ence of the quadratic curvature corrections £1 and £ 2 . 
The analytic results were known only for exact de Sit¬ 
ter inflation, and we have used the de Sitter formulas 
evaluated at horizon crossing in the case of slow-roll in¬ 


flation for which the Hubble parameter is varying. The 
validity of the method has been checked by performing 
numerical calculations. Both £1 and £2 reduce the am¬ 
plitude of primordial tensor perturbations. However, we 
have found that £2 could generate large non-Gaussianity 
of the curvature perturbation, which places a stringent 
constraint on the amount of the suppression due to £ 2 . 
Since £1 does not change the cubic interaction of the 
curvature perturbation, this evades the non-Gaussianity 
constraint. The tensor power spectrum can be as small as 
65% of the standard result due to £ 1 , which brings many 
inflation models with large tensor modes to the observa- 
tionally preferred region in the n^-r plane. We have seen 
that the tensor tilt is also modified, though the spectrum 
can never be blue. 
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